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Minkowski Functionals are a powerful tool for analyzing large scale structure, in particular if the 
distribution of matter is highly non-Gaussian, as it is in models in which cosmic strings contribute 
to structure formation. Here we apply Minkowski functionals to 21cm maps which arise if structure 
is seeded by a scaling distribution of cosmic strings embeddded in background fluctuations, and then 
test for the statistical significance of the cosmic string signals using the Fisher combined probability 
test. We find that this method allows for detection of cosmic strings with G/i > 5 x 10~*, which 
would be improvement over current limits by a factor of about 3. 



I. INTRODUCTION 

As a consequence of the realization that in many infla- 
tionary models (both supergravity-based T\\ and string- 
based PJ) the period of inflation ends with the formation 
of a network of cosmic strings, interest in searching for 
the signatures of such strings in cosmological observa- 
tions is increasing [53] . If they are topologically stable, 
the network of cosmic strings will persist at all times 
and will approach a 'scaling solution' characterized by a 
string distribution which is statistically independent of 
time if lengths are scaled to the Hubble radius (see e.g. 
[3H5] for reviews of cosmic strings and their consequences 
for cosmology). 

Cosmic strings are linear topological defects which 
arise in a wide range of quantum field theories as a con- 
sequence of a symmetry breaking phase transition. They 
are lines of trapped energy density and tension which is 
equal in magnitude to the energy density. They are anal- 
ogous to defect lines in crystals or to vortex lines in super- 
conductors and superfluids, except for the fact that the 
cosmic strings arise in relativistic field theories and hence 
obey relativistic dynamical equations of motion as op- 
posed to the strings in condensed matter systems whose 
dynamics is typically friction-dominated. The trapped 
energy density associated with the strings leads to con- 
sequences in cosmology which result in clear observa- 
tional signatures. Since the energy density per unit string 
length increases as rf' as the symmetry breaking scale 77 
increases, cosmology provides the ideal venue to search 
for strings at very high energy scales. Thus, searching 
for strings in cosmological observations provides an av- 
enue complementary to accelerator experiments for look- 
ing for signals of physics beyond the Standard Model: 
accelerator signals are more easily seen for low symme- 
try breaking scales, while signals in cosmology are more 
easily detected for larger values of 77. 

Because of the scaling solution of the cosmic string net- 
work, cosmic strings lead to a scale-invariant spectrum 
of cosmological perturbations, as in inflationary mod- 
els (see e.g. [7]). However, the string- induced density 
fluctuations are highly non-Gaussian, an effect which is 
eliminated when computing the usual power spectrum 



of density fluctuations. It follows that string signals are 
much easier to detect in position space than in momen- 
tum space. 

The scaling network contains two types of strings, a 
network of 'infinite' strings [54 and a distribution of 
string loops with radii smaller than the Hubble radius. 
String loops oscillate because of their relativistic ten- 
sion and slowly decay by emitting gravitational radiation. 
Their gravitational effects are similar to those of a point 
mass [7]. Since according to numerical simulations of 
cosmic string networks ^ the long string network dom- 
inates, we will in this paper not further consider string 
loops. 

Due to the relativistic tension of the strings, an infi- 
nite string segment will typically have a velocity in the 
plane perpendicular to the string which is of the order of 
the speed of light. Since space perpendicular to a long 
straight string is conical with a 'deficit angle' given by 

m 

a = SirGfj., (1) 

where G is Newton's constant and fi is the mass per unit 
length of the string (which is proportional to ry^) , a string 
moving with velocity Vs through the matter gas of the 
early universe at some time U will lead to a 'wake' [lOj . a 
thin wedge behind the string with twice the background 
density whose planar dimensions are 

ciU X VsjsU ■ (2) 

Here, 7s is the relativistic 7 factor associated with the 
velocity Vs, and ci is a constant of order one which is 
proportional to the curvature radius of the long string 
network in units of the Hubble radius. The mean thick- 
ness h of the wake is initially 

h = AnGfiVs-fsU (3) 

and for t > ti it increases by gravitationally accreting 
matter from above and below the wake [TT] . 

Current observations (in particular the oscillatory fea- 
tures in the angular power spectrum of cosmic microwave 
background (CMB) anisotropies) set an upper bound on 
the string tension of the order G^i < 1.7 x 10""^ [HI [13], 



2 



which corresponds to cosmic strings contributing less 
than 5% to the total spectrum of inhomogeneities. We 
will discuss this in more detail in sectionjll] In this paper, 
we are therefore considering a setup in which there is a 
contribution of cosmic strings to the power spectrum in 
addition to the dominant component of Gaussian, nearly 
scale-invariant fluctuations (such as can be produced in 
inflationary cosmology [131 [H] or in string gas cosmology 
0). 

Cosmic string wakes give rise to distinctive signatures 
in the large-scale distribution of matter in the universe. 
As discussed in [TB], they lead to planar structures in 
the distribution of galaxies. Since wakes present between 
the time of last scattering of the CMB and today rep- 
resent overdense regions of electrons, they lead to extra 
scattering of CMB photons and hence to distinctive rect- 
angular regions in the sky with extra CMB polarization 
[TT] (with statistically equal B-mode and E-mode com- 
ponents). The effect relevant to the current paper is that 
wakes represent overdense regions of neutral hydrogen 
and hence [18] lead to wedge-like regions in 21cm red- 
shift surveys of extra 21cm absorption or emission |55j . 
The amplitude of the position space signal is independent 
of the string tension and can be as large as lOOmK 
[56j . The width of the wedge, however, depends lin- 
early on fi. A crucial point is that wakes exist as non- 
linear structures at high redshifts since as soon as the 
string passes by, a wake with overdensity 2 forms. This 
is to be contrasted to the situation in Gaussian models 
with a scale-invariant spectrum of primordial cosmologi- 
cal perturbations in which no nonlinear structures exist 
until quite low redshifts. In particular, at redshifts larger 
than that of reionization, the cosmic string signal should 
stand out against the effects of Gaussian fluctuations and 
noise. Thus, 21cm redshift surveys appear to be an ideal 
window to search for cosmic string signals. 

If cosmic strings exist, then the induced 21cm maps 
will not only contain the characteristic wedges of extra 
absorption or emission. They will also contain noise, 
most importantly Gaussian noise from the primordial 
Gaussian density fluctuations which must be present in 
addition to the string-induced perturbations. Good sta- 
tistical tools will be required in order to extract the string 
signals in a quantitative and reliable way. 

In this report we study the possibility that Minkowski 
Functionals [2D] , a tool for characterizing structure which 
is orthogonal to the usual way of characterizing maps us- 
ing correlation functions, can be applied to identify cos- 
mic string signals in 21cm redshift surveys at reshifts 
larger than that of reionization. The theory behind 
Minkowski Functionals has its roots in Integral Geom- 
etry and Hadwiger's Theorem, which states that a d- 
dimensional set of convex bodies can be completely de- 
scribed by d + 1 functionals. Minkowski functionals have 
been applied numerous times in cosmology, e.g. to the 
distribution of galaxies on large scales ^21] and the analy- 
sis of CMB maps [H] . There have also been attempts to 
apply Minkowski functionals to search for signatures of 



cosmic strings in the large-scale structure of the distribu- 
tion of galaxies [23]. Since for 21cm redshift surveys the 
contributions to the 'noise' from the Gaussian source of 
primordial fluctuations is in the linear regime at the high 
redshifts we are interested in (larger than the redshift of 
reionization), we expect that Minkowski functionals will 
be very powerful in extracting cosmic string signals in 
the case of 21cm maps. 

The outline of this paper is as follows: in the following 
section we outline current observational limits on cosmic 
strings as well as prospects for future detection. In Sec- 
tion III we briefly review the toy model for the cosmic 
string scaling distribution which we use and describe the 
theoretical 21cm maps induced by cosmic strings. Sec- 
tion IV then repeats the methodology of section III but 
including a simulation of the background noise. In Sec- 
tion V we review the theory and calculation of Minkowski 
functionals, before presenting our results in Section VI. 
We conclude with a summary of our results and a discus- 
sion of future work. 



II. OBSERVATIONAL CONSTRAINTS 

The strongest constraints on cosmic strings currently 
come from measurements of the angular power spec- 
trum of the cosmic microwave background (CMB). The 
WMAP data provided no hints of cosmic strings, and 
hence placed an upper bound on the string tension [13) . 
Specifically, it constrained the string contribution to the 
primordial power spectrum to be less than 10%. How- 
ever, there are now CMB experiments with better an- 
gular resolution than WMAP provided and which can 
yield improved bounds on the cosmic string tension. The 
Atacama Cosmology Telescope (ACT) [M], a six-meter 
off-axis telescope with arcminute-scale resolution located 
in the Atacama desert in northern Chile, and the South 
Pole Telescope (SPT) [25] , a 10 meter telescope located 
in Antarctica which can study the small-scale angular 
power spectrum of the CMB, both are providing excel- 
lent data. A recent analysis of the angular power spec- 
trum of CMB anisotropics from joint SPT and WMAP7 
[T^ made use of a Markov Chain Monte Carlo likeli- 
hood analysis to place a limit on the string tension of 
G/i < 1.7 X 10~^ (at 95% confidence). This analysis was 
done in the context of the zero width cosmic string toy 
model which we will also use in the following. A similar 
limit was found by [2S], who used combined WMAP7, 
QUAD and ACT data to place limits on the tension of 
Abelian Higgs model strings. The resulting bound on the 
string tension was < 4.2 x 10^^. 

The difference in the results is mostly due to the un- 
certainties in the distribution of strings. Due to the very 
large range of scales involved (the width of a cosmic string 
is microscopic but the length is cosmological), numerical 
simulations of both field theory and Nambu-Goto strings 
involve ad hoc assumptions and/or extrapolations. The 
uncertainties in the resulting string distributions appear 
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in free parameters (e.g. the number of long string seg- 
ments per Hubble volume and the value of the constant 
Ci) of toy models distributions of cosmic strings. These 
uncertainties effect the angular power spectrum of CMB 
anisotropies. 

Local position space searches for signals of individual 
cosmic strings are less sensitive to the uncertainties in 
the string distribution than power spectra. This is the 
idea behind the proposal of [27l - f29] to search in posi- 
tion space for the line discontinuities in CMB anisotropy 
maps which long straight cosmic strings produce due to 
the Kaiser-Stebbins [SU] effect. It was proposed to ana- 
lyze position space anisotropy maps using of the Canny 
edge detection algorithm. It was shown 29 that strings 
with Gfi > 10~^ might be detectable with this method. 
Application of this method to the SPT or ACT data may 
yield interesting results. 

Another avenue for detection of cosmic strings is via 
gravitational waves, which will be produced by oscillating 
cosmic strings and/or the decay of cosmic string loops. 
Cosmic strings in fact |31| produce a scale- invariant spec- 
trum of gravitational waves with an amplitude which is 
similar to or larger than the amplitude of gravitational 
waves from simple inflationary models. Pulsar timing [32| 
or direct detection (e.g. making use of the Laser Inter- 
ferometer Gravitational- Wave Observatory (LIGO) [S3] ) 
provide means for searching for the cosmic string signal. 
However, at the moment the bounds are not competitive 
with bounds obtained from the CMB [57] . 

We are currently entering a revolution in radio as- 
tronomy, with both the Square Kilometer Array (SKA) 
[51] and the European Extremely Large Telescope (E- 
ELT)P5] aiming to be operational by early 2020. Of 
particular interest (in terms of potential to observe cos- 
mic strings) is the SKA, which as the name implies 
will consist of 1 million square meters of collecting area. 
There are many ongoing projects to develop the science 
and technology necessary to operate the SKA. These 
are categorized into: (1) Precursor Facilities, which will 
physically be at the SKA site carrying out SKA-related 
projects, (2) Pathfinder Experiments, which will develop 
SKA-related science and technology off-site, and (3) De- 
sign Studies, which will investigate technologies and de- 
velop prototypes. Many of these projects are world class 
facilities in there own right, and we will mention some of 
the experiments which may have a chance of observing 
cosmic strings. 

One such project is LOFAR [36], a Pathfinder exper- 
iment, which is measuring the neutral hydrogen frac- 
tion of the inter-galactic medium. This will map out 
the epoch of reionization between redshifts 11.5 and 6.7, 
and will be capable out of arc-minute resolution. For 
21cm measurements, the redshift range of LOFAR ex- 
tends to the "Dark Ages" before reionization, redshifts 
where the cosmic string signal will be cleanest. Another 
SKA project with potential to observe cosmic strings is 
the MEERKAT experiment [37] , which is carrying out 
a project titled MESMER that will track the neutral 



hydrogen content of the early universe by using carbon 
monoxide as a proxy. To make an accurate prediction 
as to the clumping of carbon monoxide around a cosmic 
string would require semi-analytical hydrodynamics cal- 
culations, but to a rough approximation we should expect 
to see the characteristic wedge of emission/absorption. 



III. THE COSMIC STRING SIGNAL 
A. Modelling the signal 

In any theory which leads to the formation of stable 
cosmic strings in the early universe, a network of strings 
will persist at all later times, and this network will ap- 
proach a 'scaling solution', meaning that the statistical 
properties of the network of cosmic strings are indepen- 
dent of time if all distances are scaled to the Hubble ra- 
dius. The existence of the scaling solution can be argued 
for using analytical arguments (see e.g. [IHS]), and it has 
been confirmed using extensive numerical simulations [8] . 

The scaling distribution of strings consists of a network 
of long strings of mean curvature radius C — cit (where 
Ci is a constant of order 1), and string loops with a much 
smaller radius that result from intersection and hence 
'cutting' of the long strings. As a consequence of their 
relativistic tension, the long strings have typical trans- 
lational velocities Vs close to the speed of light. Hence, 
Hubble length string segments will intersect with other 
such segments with probability close to one on a Hub- 
ble time scale. By this process, string loops are formed 
and the correlation length of the long string network 
(the mean curvature radius) increases in comoving coor- 
dianates to keep up with the Hubble radius, as confirmed 
in the numerical string evolution simulations [5]. 

As is commonly done in studies of the cosmological 
consequences of cosmic strings, we use an analytical toy 
model of the infinite string network which was first in- 
troduced in [38j : we divide the time interval of interest 
(usually the time between the time t^q of equal matter 
and radiation and the present time to) into Hubble time 
steps. In each time step we lay down an integer num- 
ber Nh per Hubble volume of straight string segments 
of length cit with random midpoints and random tan- 
gent vectors. The string segments in neighboring Hubble 
time steps are taken to be independent. The long strings 
are taken to have a velocity Vs orthogonal to their tan- 
gent vectors. Each such string will produce a wake whose 
thickness grows in time. 

We will be interested in signals which are emitted at a 
time te > ti from string wakes laid down at any time ti 
between tpq [58j and the present time . Since the planar 
dimensions of the wake (whose initial physical dimen- 
sions are given in ([2])) are constant in comoving coor- 
dinates, their physical size grows with the scale factor. 
The comoving thickness of the wake, whose initial phys- 
ical length is given by ([3]), grows linearly with the scale 
factor because of gravitational accretion. Hence, the di- 
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mensions of the wake at ig in physical coordinates are 
given by: 



^i + lfo- ^ A J. Zi + l- 

— [ciU X Uvsjs X AttGhUvsIs 



Ze + 



Z, + V 



(4) 



where is the emission redshift z{te)- 

The dominant form of the baryonic matter in the uni- 
verse before reionization is neutral hydrogen (H), which 
we detect via the 21cm hyperfine line. Measuring the 
intensity of the redshifted 21cm radiation from the sky 
has the potential of giving us three dimensional maps 
of the distribution of neutral hydrogen in the universe 
(see e.g. [311 for an in-depth review article on 21cm cos- 
mology). As explained in detail in jl8j . since wakes are 
overdense regions of neutral hydrogen, they will lead to 
excess 21cm emission or absorption. We will see this ex- 
cess in directions for which our past light cone intersects 
a wake at some time > tree, tree being the time of re- 
combination. This 21cm signal has the special geometry 
given by the wake geometry - a thin wedge in the three- 
dimensional 21cm redshift survey (see Figure 1), and it is 
this special geometry which provides the 'smoking gun' 
signal for cosmic strings, the signal which we are trying 
to extract using Minkowski functionals in this report. 

Since already at the initial time ti the wake is a non- 
linear density fluctuation, the matter being accreted onto 
the wake will collapse onto the wake, shock heat and 
thermalize [SS]. The hydrogen (H) atoms inside a wake 
laid down at redshift Zi will have a temperature at time 

given by:[T5] 



TK{te) 



167r" 



{Gnf{vsis 



, Zi + 1 m 

Zf + lkn 



(5) 



where m is the proton mass. Inserting numerical values, 
and expressing Gfj, in units of 10~^ and so written as 
(G'/i)6, we obtain 

TK{t,) « [20K]{G^i)livsJ.)'^. (6) 

Ze + l 

Note that the numerical simulations by |3T] show that the 
temperature can be taken to be approximately uniform 
inside of the wake. 

The quantity of interest is the brightness temperature 
Tf, due to 21cm transitions of the Hydrogen in the wake, 
or more specifically the difference STh in brightness tem- 
perature between photons from the wake and those from 
the surrounding space. From this point onwards, 'bright- 
ness temperature' will refer to STi,, defined by 

Tb{Ze) - T^{Ze) 



STbiZe) 



(7) 



where T^{Zf,) is the redshifted CMB temperature, and 
T&(ze) is the brightness temperature due to 21 cm emis- 
sion in the wake. A full derivation which can be found in 
[l8] yields the following expression for 6Ti,: 



Snize) = [0.07K]- 



1 



Tk 



Vl + Ze 



(8) 



where several constants have been absorbed into the pref- 
actor, which can be found in the original derivation by 
[T5] . Note that the collision coefficient Xc is given by [35] ■ 



(9) 



where T*, taken to be 0.068 K, is the temperature cor- 
responding to the hydrogen energy splitting Eiq, Aiq is 
the spontaneous emission coefficient of the 21cm transi- 
tion, n is the number density of hydrogen atoms, and 
K^Q^ is the de-excitation cross section which is given in 
[35] . Using equation ([S]), the brightness temperature for 
a given wake with a given value of z^ can be obtained, 
which in conjunction with the spatial dimensions of the 
wake, allows a 3-D map of brightness temperature to be 
generated. 



B. Generating Temperature Maps 

In the following we will outline the steps in the con- 
struction of 21cm redshift maps in the case of primordial 
perturbations produced exclusively by strings. We will 
consider hypothetical sky maps covering a patch of the 
sky of angular scale 9x9. First, we must compute the 
number of strings which contribute to the sky signal in 
this patch. This is the number of strings whose wakes 
at some point in time between tree and to intersect the 
observer's past light cone with opening angle 9. 



1. Distribution of Strings 

According to our model of the distribution of a scaling 
string network we divide the time between teq and the 
present time to into Hubble time steps. For each such 
time step centered at time ti we now compute the number 
N{ti) of string wakes which at some point in the future 
of ti will intersect the past light cone corresponding to 
an observational angle 9 (a box of size 9 x 9 on the sky) . 

We first note that N{ti) is given by the ratio of the 
comoving volume of the past light cone 



Vp{9) ^ 27(— ) ti (10) 
divided by the comoving Hubble volume at time ti 

VH{t^) = t,4, (11) 



multiplied by the number Nh of strings per Hubble vol- 
ume in the scaling solution (a number which according 
to numerical simulations is in the range 1 < Nh < 10 
fS]. This gives 



3/2 



(12) 



Now that we know the total number of string wakes 
produced at redshift Zi which will be seen in the simu- 
lation box, we need to determine the distribution of the 
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FIG. 1; The right panel depicts a sketch of the string wake-induced wedge of extra emission/absorption in 21cm redshift 
surveys. The horizontal axis represents the two angular directions, the vertical axis redshift. The left panel is a space-time 
sketch showing the position of the string which gives rise to the wedge of the right panel. Here, the vertical axis is time, and the 
horizontal axis corresponds to comoving spatial coordinates. Note, in particular, the positions where the string wedge intersects 
the past light cone of an observer at the current time tg. The almost horizontal lines represent the string at the time when it 
forms the wake (ti), and the wedge at the times si and S2 when the back (and front) of the wake intersect our past light cone. 
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redshift at which they will cross the past light cone of 
the simulation region in the sky. This can be determined 
by remembering that light travels in a straight line in 
conformal coordinates (comoving spatial coordinates and 
conformal time t). Hence, the number of wakes which in- 
tersect the past light cone at conformal time r (where the 
conformal time is measured backwards, meaning r = at 
the present time and t increases as we go back in time) 
in a conformal time interval dr is 



(13) 



where J\f is determined by demanding that the integral 
over r from Ti = T(ti) to T(io) gives Ni. The result can 
be integrated from to Te to yield the number of string 
wakes N.^^t^) which intersect the past light cone at a time 
later or equal to Te'- 



N,, 



(14) 



To determine the redshift distribution of emission times 
Ze, we need to relate the conformal time Tg to the emission 
redshift Zg. We start by expressing Tg in terms of t^,, 



d^e 

a{te) ' 



(15) 



which evaluates to 



3in 



1 - (z, + 1)- 



(16) 



An expression for the number Ni[Zf,) of string wakes with 
an emission redshift larger or equal to Zg can now be ob- 
tained by substituting the preceding equation into equa- 
tion |14[ and using our expression for Ni from |12[ This 
gives 



iV,(zg) = NH27{-^fz{Uf'^ 



1- (^e + 1)- 
1 - (z, + 1)- 



(17) 



Due to the discrete nature of this simulation, we need 
to approximate this relationship by splitting the redshift 
space into slices, and laying down strings in each slice. 
The smallest possible slice would be the number of dis- 
crete values of redshift in the simulation (128, as we will 
discuss shortly) divided by the range of redshift, which 
we shall take to be 50 for this example. However, as we 
must have an integer numbers of wakes, a slice this small 
would result in the majority of slices having wakes. 
Hence we instead choose a slice thickness Azg = 1, which 
will contain a number of wakes given by 

A7V,(Zg,Ze-l) = [1.41 X 10^] 



(18) 



-3/2 
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FIG. 2: Distribution of the redshifts at which the string wakes 
intersect the past light cone (for a fixed time Zi = 10'' at which 
the strings were produced). The horizontal axis is the emis- 
sion redshift Ze, the vertical axis is the number ANi{ze, Ze — 1) 
of string wakes produced at redshift Zi (taken to be 10"^) which 
intersect the past light cone of the sky area we are considering 
in the redshift interval between z^ — 1 and z^. 
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where we have included all prefactors into the constant. 
A plot of this distribution is shown in Figure 2. As is 
apparent, the largest number of wakes intersect the past 
light cone at low redshifts. This can be understood since 
the low redshift range covers most of the physical volume. 



2. Angular and Redshift Scales of a Wake-Induced Wedge 

As discussed in Section 2, each string wake which inter- 
sects the past light cone of the observer's angular patch 
will lead to a thin wedge of extra 21cm absorption or 
emission. We are concerned with string wakes produced 
at time ti which intersect the past light cone of the obser- 
vational patch at some time te > tree- Each such string 
wake will lead to a thin wedge in a 21cm redshift survey 
of a certain angular extent and certain thickness. 

First, let us consider the angular extent of the wedge. 
A string laid down at a time ti will create a wake with 
physical dimensions given by equation[4] To calculate the 
angular size of a wake, one must start with the physical 
length of the wake at time U which in direction tangent 
to the string is: 

XpiU) = ciU , (19) 

where (as we recall from before) ci is a constant of or- 
der 1. In the direction of string motion the size is the 
same except that the factor ci gets replaced by Vgjs- In 
our work we choose the rather realistic parameter values 
where the two constants are the same. The correspond- 
ing comoving distance is then given by 

2 

Xcik) = (^j;^ ' Xp{t,) = {z, + ly^^^Cito . (20) 

Hence, the corresponding angle 6i is: 

e, = {z, + l)-^''^ci X 90°. (21) 



For the most numerous and thickest wakes, Zi + \ = 1000, 
and if the constant ci is taken to be 1/3 , the answer is 
then approximated by: 

e, = {z, + l)-i/2 X 30° (22) 

which yields approximately 1° as the angle of the wake- 
induced wedge. 

Since the string which produces the wake is moving, 
the string-induced wake is slightly "tilted" in the three- 
dimensional redshift-angle space, i.e. the tip of the wedge 
(which corresponds to the string at the latest time) is at 
a slightly larger redshift than the mean redshift of the tail 
of the wedge (where the string was at the initial time). 
This does not change the fact that the projection of the 
wedge onto the angular plane has the dimension given 
above. However, when computing the projection into the 
redshift direction we have to be careful. The wedge is thin 
in direction perpendicular to the plane of the wedge, but 
this perpendicular direction is at an angle relative to the 
redshift direction. In the following we will compute the 
thickness of the wedge in perpendicular direction to the 
wedge plane. 

We now calculate the thickness of the wake in the per- 
pendicular direction. The finite thickness in redshift di- 
rection originates from the fact that photons originating 
from different parts of the wake are redshifted by slightly 
different amounts. The starting point of the computation 
is the formula 

valid in the matter-dominated period. Taking differen- 
tials, and setting t — yields 

Az= (^^)v'/\-V3Ai,, (24) 

where Ate is the time delay between photons from the 
top and the bottom of the wedge. We compute At^, at the 
midpoint of the wake (point of medium thickness) . This 
is given by the thickness of the wake, which, as discussed 
in il8i grows linearly in comoving coordinates because of 
gravitational accretion onto the initial overdense region. 

Using linear perturbation theory, the width of the wake 
at time tg is 

wite) = A7^G^iv,J,to{z{t,) + l)'^'(z(t2) + 1)"' , (25) 
which equals Atg. Hence, from ( p4| ) it follows that 

Az = ^G^ivsj^zit,) + l)'/'(z(<2) + 1)"'^' . (26) 

The analysis of gravitational accretion performed using 
the Zel'dovich approximation instead of with naive linear 
cosmological perturbation theory yields the same result 
[l8] except that the coefficient is 247r/5 instead of Sir/ 3. 

Once the dimensions of the wake are calculated, it re- 
mains to map this onto a 3-D lattice of points from which 
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the Minkowski Functionals can be calculated. To do this 
requires setting the resolution with which the wakes will 
be studied, and then scaling up the dimensions of the 
wake. For the purposes of this simulator, we use a lat- 
tice of side length 128, and study a volume that subtends 
an angle of 10° in both the angular directions, with Zg 
ranging from 5 to 50, and — 'i. In this case, the 
angular resolution is 10°/128 ~ .078°, while the resolu- 
tion in the redshift direction is 45/128 = .35. To scale 
up the dimensions of the wake only requires multiply- 
ing each dimension by the inverse of the resolution (ie: 
128/45 = 2.84 in the direction). 



3. Brightness Temperature: STi, 

The final step in the calculation is to determine the 
brightness temperature at every point of the three dimen- 
sional angle-redshift map {9,(j),z). Any point for which 
the past light ray in angular direction {9, (j)) does not 
intersect a wake at redshift z yields zero brightness tem- 
perature. A point which at redshift Zg is intersecting a 
wake due to a string present at time ti is assigned a non- 
vanishing brightness temperature given by Equation ^ , 
with Tpc obtained from Equation Overlapping wakes 
are assumed to be non-interacting. 



Instead, we focus on the contribution of the inhomo- 
geneities mentioned above which lead to the formation 
of mini-halos. This contribution has been studied done 
both semi-analytically, [32]: and later using large scale 
numerical simulations |43| . 

To calculate the signal from minihalos we must find the 
mass function dn/dM, which can, when integrated over 
a range of masses, give the number density for minihalos 
in this mass range. There are two possible mass profiles 
for minihalos, the Sheth-Tormen model and the Press- 
Schechter model. Recent numerical simulations support 
the Press-Schecter model, (see Figure 2 of [13]). Once 
we know the mass function, the brightness temperature 
is then calculated using line of sight integrals given a 
realization of the inhomogeneities described by the mass 
function. 

The Press-Schecter function for the number density of 
minihalos is given by |44) : 



N{M)dM = - 



1/2 



Sc da 
~^dM^ 



'dM (27) 



where p is the mean (baryonic and dark) matter density 
of the universe, and a'^{M,t) — {5\j){t) is the root mean 
square mass fluctuation on the mass scale M. Halos will 
form in regions with a 5m above a critical fluctuation 
5c « 1.69. Our number density simplifies to: 



IV. BACKGROUND NOISE 
A. Modelling the Signal 

As discussed in the introductory section, in our setup 
the cosmic strings only make a contribution of less than 
5% to the total power spectrum of inhomogeneities. The 
dominant contribution is in the form of an approximately 
scale-invariant spectrum of Gaussian perturbations such 
as can be produced in a number of cosmological scenar- 
ios, e.g. in inflationary cosmology [15] or in string gas 
cosmology [3]. The dominant Gaussian fluctuations are 
in the linear regime until late times on large length scales 
and are hence not expected to contribute a lot to 21cm 
fluctuations at high redshifts (redshifts before reioniza- 
tion). There will, however, be smaller scale fluctuations 
which become non-linear, forming so-called "mini-halos" 
which then contribute to the spectrum of 21cm fluctua- 
tions. It is crucial for us to check that the cosmic string 
effects can be detected above the noise from the Gaus- 
sian fluctuations (which we call "background noise" in 
the following). 

The contribution of "background noise" to the 21cm 
fluctuations of the universe has been studied in detail 
(see e.g. [121133]). As was shown, there is an effect on 
21cm maps which comes from the diffuse inter-galactic 
medium (IGM). However, this effect is homogeneous in 
space on the cosmological scales which interest us here 
and will hence not be further considered in this paper. 



N{M,z)dM = a, 



(3-f n)\ / p 



6 



2V'^ f M 



X exp 



Try \M^{z) 

1 / M ^ (3+n)/3' 



2 \M^{z) 



(28) 
dM 



where is the mass scale for which the fractional den- 
sity fluctuation equals 5c, and where we assume a primor- 
dial power spectrum of density fluctuations P{k) ss fc" 
(here, P{k) does not contain the phase space factor 
and is not dimensionless. A scale-invariant spectrum of 
curvature fluctuations corresponds to n = 1). 

The spatially averaged brightness temperature of a set 
of halos is given by |:42! : 



5% 



0(1 + zY 
vqH{z') Jm, 



dN 

l^Vcfs5n..,A—dM (29) 



where ^Veff is the redshifted effective line width, ST},^,^^ 
is the brightness temperature of individual minihalo at 
the frequency i/q, and A is the geometric cross section 
of the minihalo. The integral (worked out in (35] - see 
Figure 3 in that paper) yields a noise temperature which 
peaks at a value of 4mK at a redshift of 10 and decays 
at higher redshifts. 
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B. Generating Noise Temperature Maps 

Given that the signal from the IGM is spatially uni- 
form, we focus on the contribution from minihalos. 
The corresponding noise temperature maps are con- 
structed by modeling the spatial fluctuations as a three- 
dimensional Gaussian random field with a power spec- 
trum Prik) (again defined without the phase space factor 
k^) with a slope which will be discussed below, and with 
a variance given by ( 29 1 . The amplitude of the resulting 



noise is rescaled in redshift direction by incorporating the 
redshift dependence of the variance from ( 29 ) . 



To accomplish this, we follow the procedure laid out 
in |29j . which we generalize to 3 spatial dimensions. We 
can expand the spatial fluctuations in temperature into 
hyperspherical harmonics.: 



AT 



= X! V'), (30) 

Lm.n 



where the Yim,n are the spherical harmonics generalized 
to 3-t-l dimensions, the ^hyperspherical harmonics\ and 
the ai_rn,n are the coefficients of this expansion. We can 
simplify this expression to a sum of plane waves, by mak- 
ing the flat sky approximation (valid for angular scales 
< 60°). This gives us: 



fix) = ^f(^)e 



(31) 



where we have used the abbreviation T{x) 
Comparing this decomposition with our previous expres- 
sion in terms of hyperspherical harmonics, we see that 
the T{k) correspond to the a/^m,n and hence the T{k) 
give us the power spectrum: 



Prik) 



(32) 



where the temperature power spectrum Pt is defined by 
the last equality. 

Since the mini-halos are produced by the density fluc- 
tuations we use the power spectrum of density fluctua- 
tions to determine the slope of Prik). Thus, we take 
Pxik) to be proportional to the primordial power spec- 
trum: 



Prik) oc fc" with n = 1 , 



(33) 



and the amplitude is determined by demanding that the 
variance is given by (29). First, in fact, we generate 



a three-dimensional Gaussian random field Tqrf with 
variance 1. 

Note that in a model with a primordial scale-invariant 
spectrum of curvature perturbations, the dimensionless 
power spectrum of fractional density perturbations is ap- 
proximately scale-invariant on scales which entered the 
Hubble radius before the time of equal matter and ra- 
diation. This is a consequence of the processing of the 
primordial spectrum which happens on sub-Hubble scales 



(see e.g. any text which discusses the Newtonian theory 
of cosmological perturbations, e.g. |46]). The spectrum 
on small scales thus takes on a slope n = —3. We do not 
consider this effect in the current simulations. 

Using the ergodic hypothesis, we construct a Gaussian 
noise field by drawing the T(k) for fixed k = \k\ from 
a Gaussian distribution with power given by Pxik) and 
amplitude chosen to give variance 1. 

We do our calculations on a lattice of (fci, A:2, fca) coor- 
dinate values ranging from to N^ax — 1- We convert to 
the corresponding k values with: 



277, ~ 



L 



(34) 



where k^ax corresponds to the angular resolution of the 
lattice in the ki direction. 

We can then calculate the T{k) using: 



f{k) 



^|^^(g,{k)+^g,(k)) (35) 



where gi and g2 are randomly generated from a Gaussian 
distribution with variance 1 and mean 0, and we enforce 
f{k) = f{-k) to ensure that T{x) is real. 

We can then construct a spatial map by taking the 
inverse Fourier transform of the preceeding expression. 
The result of this is a Gaussian random field Tqrf with 
a correlation function specified by the power spectrum 
and variance 1. We now wish to enforce the redshift 
dependence of the amplitude, which we do by identifying 
redshift z with the third spatial direction. Hence the final 
noise field is specified by: 



Tnoise{x,y,z) = 5Ti,{z)Tgrt{x , z) 



(36) 



with amplitude 5Th{z) given by (29). Note that since the 
change in amplitude is uniform over the angular direc- 
tions, this introduces deviations from strict Gaussianity 
of the the three-dimensional distribution. 



V. INTEGRAL GEOMETRY AND MINKOWSKI 
FUNCTIONALS 

Minkowski Functionals are a useful tool to analyze the 
topology of a d-dimensional map. Originally, these func- 
tionals were considered in the context of describing the 
topology of a body B embedded in a d-dimensional Eu- 
clidean space using an approximation in which the body 
B was approximated by the set of convex bodies of a 
particular form but varying size (hence we have func- 
tionals and not just functions). In our application we will 
consider the Minkowski functionals as functionals of the 
class of bodies (termed an 'excursion set') which enclose 
regions of the map where the value of the variable which 
is mapped is larger than a cutoff value and we consider 
varying the cutoff value. For example, applied to CMB 
temperature anisotropy maps we consider bodies which 
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are regions of the sky where the temperature T is larger 
than a cutoff temperature Tc whose amphtude we vary. 
In the case of interest here we consider three dimensional 
21cm brightness temperature maps and we consider the 
topology of the volumes which contain points where the 
temperature exceeds a critical temperature whose mag- 
nitude we vary. 

Minkowski first developed these functionals ^47, in 
1903 to solve problems of stochastic geometry [JS]. This 
then led to the development of Integral Geometry in the 
mid 1900's. At the heart of Integral Geometry is Had- 
wiger's Theorem |2Q], which deals with the problem of 
characterizing the topology of the body B using scalar 
functionals V . These functionals must satisfy certain re- 
quirements [49] : 

1. Motion Invar iance : The functionals should be 
independent of the position and orientation in space 
of the body. 

2. Additivity: The functionals applied to the union 
of two bodies equal the sum of their functionals 
minus the functionals of their intersection: 

V{A UB) = V{A) + V{B) - V{A n B) . (37) 

3. Conditional Continuity: The functionals of con- 
vex approximations to a convex body converge to 
the functionals of the body. 

Hadwiger's Theorem [20J states that for any d dimen- 
sional convex body, there exist d -I- 1 functionals that sat- 
isfy these requirements, denoted V},j = 0,...d, for the 
j'th functional. Furthermore, these functionals provide 
a complete description of the body. Mathematically, the 
j'th functional of a d-dimensional body B is an integral 
over a (d — j) -dimensional surface of B. For example, 
in three dimensions, the first functional Vb is simply the 
volume of the body, and the second (Vi) is the surface 
area. In all dimensions, the last functional Vj is given by 
the Euler characteristic x, which in three dimensions is 
defined gH] as: 

X — number of components - number of tunnels 

+ number of cavities . (38) 

The simplest example is a set of balls of radius . When 
Tj is very small, there are no intersections of balls, and 
X is very close to the number of balls. As increases, 
the balls will begin to intersect, creating tunnels. Thus 
X will become negative. At a certain threshold value of 
r^, X "will once again be positive, as tunnels arc cut off to 
form cavities. Finally, when becomes very large, the 
entire space is filled, and x is equal to 0. 



The geometric meaning of the functionals for d < 3 is 
summarized in Table I [35] . A key feature of Minkowski 



TABLE I: Geometric Interpretation of Minkowski Functionals 
in 1,2, and 3 dimensions [IS] 



d 1 


2 3 


Vb length 


area volume 


Vi X 


circumference surface area 


V2 - 


X total mean curvature 


14 - 


X 



Functionals is that they incorporate information about 
correlation functions of arbitrary order ^48j . Hence, the 
Minkowski Functionals are much more sensitive to non- 
Gaussianity then a three or four point correlation func- 
tion, making them an ideal tool to search for signatures of 
topological defects such as cosmic strings. An important 
fact is that exact expressions exist for a Gaussian Field 
(given by Tomita in [50] ). allowing for a visual compar- 
ison between string-based and Gaussian maps (since all 
we have to look for is statistically significant deviations 
from the curves for a Gaussian distribution). We will 
apply this method in Section 6. 

Next we will apply Minkowski functionals to the tem- 
perature maps whose construction is described in the pre- 
vious section. For any fixed value of the temperature, we 
consider the corresponding iso-temperature surface. The 
Minkowski functionals probe the topology of this set of 
surfaces (or equivalently the bodies they enclose). They 
are calculated as a function of the temperature threshold, 
and the results are plotted as a function of this threshold. 
It is standard practice to use as the x axis not the actual 
value of the field, but a variable v — f{x)/cr which is the 
number of standard deviations a by which the value of 
the field f{x) (scaled to have zero spatial average) de- 
viates from zero . Hence, a rescaled threshold v = ±1 
implies a fluctuation of the field by one standard devia- 
tion from the average value. The reason for using v is 
to remove the effect of a constant factor on the function- 
als, and hence to allow for a comparison with structure 
that may be identical in every way but the magnitude of 
the temperature. For the purposes of this simulator, we 
calculated the functionals for the range —4 <v<A. 

The calculation of the Minkowski functionals was 
done using the program MinkowskiS, written by Thomas 
Buchert. This piece of software takes a 3-D array of fioats 
as input (in binary format), and outputs two estimates 
(with error values) of the Minkowski Functionals as well 
as the functionals for a Gaussian field. The two estimates 
for the Minkowski functionals correspond to two meth- 
ods, one derived from differential geometry and the other 
from integral geometry. As a full derivation of each can 
be found in [51 , only a brief explanation will be provided 
here. 

In the context of differential geometry, it is possible to 
describe all local curvature in terms of geometric invari- 
ants (Koenderink Invariants) . We then express the global 
Minkowski functionals Vk in terms of the local Minkowski 
functionals V^'"*^-* , which are calculated in [5T]. For the 
three dimensional case, the functionals of a field i^(x) in 
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a volume T) are given by 



B. Functionals of Noise 



1 

W\ 



d'x5{u - v{-k)) |Vi.(x)| Vj:^°''\v,^) . 



(39) 

In the discrete case, this is equivalent to summing over 
the lattice and taking a spatial average. 

In the context of integral geometry, the Minkowski 
Functionals are calculated using Crofton's Formula. This 
provides a very elegant expression for the fc'th functional 
of a d dimensional object K, in a volume that consists of 
L lattice points of cubic lattice spacing a: 



1 



fc 



d-k+j 



(40) 

where u!j is the volume of a j-dimensional unit sphere, 
and Nj{K) is the number of j-dimensional lattice cells 
contained in K. For example, Nq{K) is the number of 
lattice points in K, and N3{K) gives the number of cubes 
contained within K. 



VI. RESULTS 
A. Functionals of Cosmic Strings 

We first consider the Minkowski Functionals of the 
wake signal alone, by taking a scaling solution of strings 
with (G/i)6 = 0.17, using the notation (G/i)6 = G/i x 
10^^. The results are shown in Figure 3. The four boxes 
show the mean Minkowski functionals after 100 simula- 
tions, with the errors taken to be the standard deviation. 
In each box, the green curve is based on calculating the 
functionals using the method of Koenderink Invariants, 
and the red curve is based on using the Crofton Formula 
(the red and green curves are almost perfect replicas of 
one another, and so the green is only visible in regions 
where they differ). 

Since the number density of wakes is peaked around 
high redshift, the mean brightness temperature fiuctua- 
tion 6Tf, will be negative, and hence the threshold value 
1/ — corresponds to a slightly negative brightness tem- 
perature fluctuation. The threshold value Vc correspond- 
ing to STh = is a small positive number. It is at 
this threshold value that the four Minkowski function- 
als change abruptly. For v > Vc the volume with STh 
greater than the threshold value vanishes, and hence the 
functionals Vi , V2 and V3 vanish. For v < i^c the high 
temperature volume does not vanish, but it corresponds 
to the outside of the wakes as opposed to the inside. 
Hence, the integrated mean curvature is negative. The 
key feature to notice in these plots is the abrupt change 
that occurs at the threshold value and the asymmetry 
about this threshold. 



We now turn our attention to the Minkowski Function- 
als of the background noise due to minihalos, shown in 
Figure 4, where we again calculate the average function- 
als over 100 simulations and indicate the standard devi- 
ation with error bars. The corresponding functionals for 
a Gaussian random field are shown as the dashed black 
curves, whereas the functionals for the noise maps are 
shown in red and green. The key feature of these plots is 
the symmetry (or antisymmetry in the case of V2 and Vq) 
about 0, which matches the results for a Gaussian ran- 
dom field. However the noise map functionals all share 
the common characteristic of being tightly concentrated 
,about 0, which allows the noise map to be differentiated 
from a pure Gaussian random field. As mentioned above, 
the deviation of the noise maps from being a pure Gaus- 
sian random field comes from the scaling of the ampli- 
tude with redshift which is uniform over the two angular 
coordinates. 



C. Functionals of Cosmic Strings Embedded in 
Noise 

Armed with our knowledge of the behaviour of the 
Minkowski Functionals for both pure cosmic string wake 
and background noise maps, we can start with the real 
fun: differentiating maps of wakes and noise from maps 
of pure noise. To do this, we generate a map of wakes 
and noise by adding (at every point in space) the tem- 
perature fiuctuations from wakes and the fiuctuations 
from the mini-halo noise maps. We can then compare 
the Minkowski functionals as we vary the string tension, 
allowing us to place constraints on the minimum string 
tension necessary for the strings to be detected via this 
method. 

Before moving on to plots, let us develop our method 
for testing significance. Consider a threshold bin v = vq, 
which corresponds to a mean wake [60j functional and 
mean noise functional Vn, with standard deviations 
and cr„. This allows us to define a new variable d which 
is the difference between the sets: 



d^Vy,-Vn 

such that: 

d = V^-Vn 
with a standard error given by: 



(41) 
(42) 

(43) 



where n is the number of simulations. The corresponding 
t — value is defined by: 



t = 



SEd 



(44) 
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FIG. 3: Minkowski functionals of a pure scaling solution of cosmic string wakes, with G/j, = 0.17 x 10~ . Each of the four 
boxes shows the results for one of the four functionals, averaged over 100 simulations each using a lattice of 128"^ pixels. Each 
box contains two curves, corresponding to computing the functionals computed making use of one of the two methods, either 
Koenderink Invariants (the green curve) or the Crofton Formula (the red curve). 



1.5 





Since t-values follow a t-distribution, we can convert this 
into a, p — value : 



p = 2{l-tcdf{t,df)) 



(45) 



into a using the Fisher combined probability test: 



(47) 



where tcdf is a cumulative density function of a t distri- 
bution, the factor of 2 is because we want the 2-tailed 
distribution, and df is the number of degrees of freedom, 
given by: 



df = {n - 1) 



(46) 



where n is the number of data points in each of the sets 
and Vn (the number of maps that contributed to the 
mean values). By applying this at every threshold bin, 
we can calculate a set of p- values, which can be combined 



where k is the number of p- values being combined. Note 
that this will follow a distribution with k degrees of 
freedom. 

The important thing to note in this procedure is that 
the is sensitive to the standard error, as opposed to 
solely the standard deviation of the functionals. Hence, 
a more fitting visual indicator of the error on the func- 
tionals is a I \fn. 

A comparison of the functionals of (wakes-|-noise) 
against pure noise is shown in Figures 5 to 8, which 
were done using a lattice size of 128^, each for a dif- 
ferent value of G/x. In each set of graphs, the black 
curves represent the results for the pure noise maps, and 
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FIG. 4: Minkowski functionals of tlie background noise, averaged over 100 simulations each using a lattice of 128 pixels. 
Each box contains three curves. Two of them are for the noise map with the functionals computed making use of one of two 
estimates, either Koenderink Invariants (the green curve) or the Crofton Formula (the red curve). The dashed black curves 
give the Minkowski functionals for a Gaussian map. 
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orange curves (the less symmetric - in the case of Vi 
and V3 - or less antisymmetric - in the case of Vq and 
V2 - curves) represent what is obtained for maps con- 
taining both strings and noise. The found to 
be infinite for G/i > 6.7 x 10^^, but converged to for 
G/i = 4.2 X 10~®. Hence from this preliminary analysis 
we conclude that the minimum string tension G/z must 
be roughly 5 x 10~^ for the Minkowski functional method 
to be able to detect the cosmic string signals which are 
embedded in the noise. 



VII. SUMMARY AND OUTLOOK 

This purpose of this investigation has been to both 
demonstrate the power of Minkowski Functionals for de- 



tecting non-Gaussian behaviour, and to investigate the 
topology of 21cm distributions in cosmological models 
with a scaling distribution of cosmic string wakes. We 
generated 3-D maps of the 21cm brightness temperature 
in pure cosmic string models, in models with only back- 
ground noise from mini-halos, and in maps in which both 
sources of 21cm fluctuations were present, using meth- 
ods similar to those employed for constructing 2-D CMB 
temperature by |29) . These maps were then be inves- 
tigated with Minkowski functionals, using the software 
'MinkowskiS' made available by Buchert [S2]. For large 
values of the string tension, the non-Gaussian signatures 
of string wakes are clearly visible in the Minkowski func- 
tionals. We studied for which range of values of the string 
tension the difference between the strings -I- noise versus 
the pure noise maps as seen in the Minkowski function- 
als are statictically significant (making use of the Fisher 
combined probability test). A large difference was found 
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FIG. 5; Minkowski functionals of a combined map (wakes+noise) with G/i = 8.5 x 10"** along with the corresponding functionals 
for a pure noise map, averaged over 100 simulations each using a lattice of 128"^ pixels. The black curves are for the pure noise 
maps, the orange ones for the strings+noise maps. 
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between the functionals of the string-induced brightness 
temperature map and the corresponding functionals of 
background noise, v^rith maps of strings embedded in 
noise found to be differentiable from maps of pure noise 
when Gfj.> 5x 10"*. 

The analysis presented in this report serves as a proof 
of concept that the Minkowski functionals can be used as 
a statistical tool to search for cosmic string signals. To 
extend this preliminary analysis would mean perform- 
ing large scale numerical simulations, as new features 
may present themselves as we study the 21cm maps with 
greater resolution. 
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There are also non-inflationary cosmological models like 
string gas cosmology in which a network of cosmic 
strings may be produced. 

The network of 'infinite' strings can be viewed as a ran- 
dom walk of strings with a step length comparable to the 
Hubble scale. 

Whether the signal is in emission or absorption depends 
on the ratio of the gas temperature in the wake and the 
CMB background temperature. 

The Fourier space amplitude, on the other hand, does 
depend on /i [19]. 

The gravitational wave signals from cosmic string cusps 
might be larger, but their amplitude is very uncertain. 
Wakes formed at earlier times have smaller angular ex- 
tent and thickness. The thickness cannot grow until t^q. 
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Hence, such wakes will be less visible than the ones we 
study. 

[59] This is true if the kinetic temperature which the hydrogen 
atoms collapsing onto the wakes acquire is larger than 
2.5Tj, where Tg is the temperature of the gas outside of 
a wake |18) . This will be true for sufficiently large values 
of Gfi. For smaller values, the wake is diffuse [40]. The 
density contrast inside the wake will be smaller, but the 



thickness will be larger. These effects will compensate 
eachother in a Minkowski functional analysis of the wake 
distribution provided that the smoothing length in the 
Minkowski code is sufficiently large. 
[60] From now on "wake" simulation will denote the wake 
plus noise simulation. 



